AMENDMENTS 
In the Specification 

[0031] This particular problem motivated us to introduce the radial frame multiresolution analysis. 
Our construction is based on a very general multiresolution scheme of abstract Hilbert spaces 
developed by Papadakis in ft6J M. Papadakis. Generalized Frame Multiresolution Analysis of 
Abstract Hilbert Spaces. 200 1 , namely the Generalized Frame Multiresolution Analysis (GFMRA). 
The main characteristic of GFMRAs is that they can be generated by redundant sets of frame scaling 
functions, hi fact, GFMRAs encompass all classical MRAs in one and multidimensions as well as 
the FMRAs of Benedetto and Li (f6^ see J J. Benedetto and S, Li. The Theory of Multiresolution 
Analysis Frames and Applications to Filter Banks. Appl. Comp. Harm. Anal, 5:389-427. 1998) . 
[0032] hi this invention, we construct non-separable Shannon-like GFMRAs of L^(E") whose scaling 
functions are radial and are defined with respect to certain unitary systems, which we will later 
introduce. We also derive certain of their associated frame multi wavelet sets. Our construction is the 
first of its kind. Scaling functions that are radial have not been constructed in the past. However, 
certain classes of non separable scaling functions in two dimensions, with some continuity properties 
with respect to dyadic dilations or dilations induced by the Quincunx matrix only have been 
constructed in the past (e.g., [7, 12, 11 8], [4] A. Cohen and L Daubechies. Nonseparable 
Bidimensional Wavelet Bases. Revista Matematica Iberoamericana. 9:51-137, 1 993; J. Kovacevic 
and M. Vetterli. Nonseparable Multidimensional Perfect Reconstruction Filter-banks. IEEE 
Transactions on Information Theory. 38:533-555. 1992: W. He and M.J. Lai. Examples of Bivariate 
Nonseparable Compactly Supported Orthonormal Continuous Wavelets. In M. Unser, A. Aldroubi, 
A. Laine editor. Wavelet Applications in Simal and Imase Processing IV, volume 3169 of 
Proceedinss SPIE. pages 303-3 14. 1997: K. Grochenig and W. Madvch. Multiresolution Analysis. 
Haar Bases and Sdf-SimilecrTiVmQsJEEE Transactions on Information Theory, 38:558-568, 1992: 
and A. Ayache. E. Belogay, and Y. Wang. Orthogonal Lifting: Constructing New (Symmetric) 
Orthogonal Scaling Functions. 2002.) . All of them have no axial symmetries and are not smooth, 
except those constructed in f5J E. Belogay and Y. Wang. Arbitrarily Smooth Orthogonal 
Nonseparable Wavelets in r^. SIAM Journal of Mathematical Analysis. 30:678-697. 1999. which 
can be made arbitrarily smooth, but are highly asymmetric. Another construction in the spirit of 
digital filter design, but not directly related to wavelets can be found in ftj E.H. Adelson. E. 
Simoncelli. and R. Hingoranp. Orthogonal Pyramid Transforms for Image Coding, hi Visual 



Communications and Image Processing IL Volume 845 of Proceedings SPIE. pages 50-58. 1987 
and fW E.P. Simoncelli, W.T. Freeman, E.H. Adelson, and J. P. Hager. Shiftable Multi-Scale 
Transforms. IEEE Transactions Information Theory. 38(2):587-607, 1992 . The latter construction 
and this of curvelets {e.g., see J. Starck, E. J. Candes, and D.L. Donoho. The Curvelet 
Transform for Image Denoising. IEEE Transactions Imase Processins:. 1 1(6): 670-684, 2002) share 
two properties of our Radial GFMRAs: the separability of the designed filters with respect to polar 
coordinates and the redundancy of the induced representations. However, our constmction in contrast 
to those due to Simoncelli et. al. and to Starck et al. are in the spirit of classical multiresolution 
analysis and can be carried out to any number of dimensions and with respect to a variety of dilation 
matrices. 

[0033] The merit of non separable wavelets and scaling functions is that the resulting processing of 
images is more compatible with that of human or mammalian vision, because mammals do not 
process images vertically and horizontally as separable filter banks resulting from separable 
multiresolution analyses do (f20J M. Vetterli and J. Kovacevic. Wavelets and Subband Coding. 
Prentice Hall PTR, Englewood Cliffs, NJ, 1995) . As Marr suggests in his book D. Marr. Vision, 
A Computational Investigation into the Human Representation and Processing of Visual 
Information. W. H. Freeman and Co., New York, N.Y., 1982 , the human visual system critically 
depends on edge detection. In order to model this detection, Marr and Hildreth used the Laplacian 
operator, which is a "lowest order isotropic operator" (B4^ D. Marr and E. Hildreth. The Theory of 
Edge Detection. Proc, /?. Soc, London B, 207:187-217, 1980) , because our visual system is 
orientation insensitive to edge detection. Thus, the most desirable property in filter design for image 
processing is the isotropy of the filter. Radial scaling functions for multiresolutions based on frames 
are the best (and, according to proposition 5, the only) type of image processing filters that meet the 
isotropy requirement. 

(0041) There exists a mapping uiG-^ G satisfying 

gD = /)a(g), for every geG 
This particular assumption implies that a is an injective homomorphism and a(G) is a 
subgroup of G, (See f95 D. Han, D.R. Larson, M. Papadakis, and T. Stavropoulos. 
Multiresolution Analysis of Abstract Hilbert Spaces and Wandering Subspaces. In D. R. 
Larson L. Baggett, editor. The Functional and Harmonic Analysis of Wavelets and Frames, 



volume 247 of Cont. Math., pages 259-284. Amer. Math. Soc. 1999 for proofs) 
\G : a(G)\ = n< +«>, where \G : a{G)\ is the index of the subgroup a(G). 

[0062] The frame scaling function can be determined in terms of Bessel functions, because it is a 

radial function. 

The proof of equation (7) can be found in \¥fj M. Pinsky. An Introduction to Fourier Analysis and 
Wavelets, 2001 Lemma 2.5.1 . 

[0063] We will not give any details regarding Bessel fiinctions. However, the reader may refer to 
M. Pinsky. An Introduction to Fourier Analysis and Wavelets, 2001 and G.E. Andrews, R. 
Askev, and R Roy. Special Functions, Number 71 in Encyclopedia of Mathematics. Cambridge 
University Press, 2000 for an extensive treatment of their main properties and of course to the bible 
of the topic fii] G.N. Watson. A Treatise on the Theory of Bessel Functions. Cambridge 
Mathematical Library. Cambridge University Press, 1944. H ere, we only include the following 
formula. 

The function, J„, given by the above equation is called the Bessel function of the first kind of order 
a. 

(0078] Using a theorem relating to the radial FMRA {Vj}j, we first set f^y = Fiy^) 

Wj-FiWj) ,wherey6if^o = H^o) = ^^^^ 

a unitary operator on £^(1"). Combining these facts with = I^(B^^( D)) 

where Qis the annulus D fl (^'(D))*", and the superscript c denotes the set-theoretic complement. 
Since an arbitrary orthogonal projection R defined on a Hilbert space H maps every orthonormal 
basis of H onto a Parseval frame for R(H) { [2, 10] A. Aldroubi. Portraits of Frames. Proceeding of 
the American Mathematical Society, 123: 1661-1668, 1995, and D. Han and D.R. Larson. Frames. 
Bases and Group Representations, volume 147 o{ Memoirs, American Mathematical Society, 2000) , 



we obtain that the orthogonal projection defined on L^{T") by multiplication with the indicator 

function of Ogives a Parseval frame for Z,^(Q, namely the set {e^Xo : k 6 Z"}. 

[0082] One of the instrumental tools of this construction is the square root of the autocorrelation 

function <D, which is defined by^(^y = O(^), a.e. on T", Also, the inverse of/i(Q is defined on the 

range of O(^) and is denoted by It can also been proved that the range projection P of the 

Analysis operator 5 is defined by Pco(^) = P(^)a)(^), where o) e /.^(T"), and that for a.e ^ 6 the 

range projection of <1>(^) is the projection P{1Q, For the sake of completeness, it must be noted that 

P{') is a projection- valued weakly measurable function defined on T^. Since <i> = , we deduce 

= Xoi^ T^ The latter observation in conjunction with the preceding argument 

implies that A{1Q'^ = 1 , if ^ 6 D. For all other ^ e T^ we have A{i^) = 0, so for these ^, we adopt the 
notational convention = 0. Last but not least, an abelian group very instrumental in the 

discussion that follows is the kernel of the homomorphism p defined by 

The latter equation implies that, for every 6 T^ p(^) is the unique point in T^ such that p(^) + k 
= A'^t^ . The kernel of p is homeomorphically isomorphic to a dual group of the quotient group 
DiAiJ}) as shown in greater detail in reference M. Papadakis. Generalized Frame 
Multiresolution Analvsis of Abstract Hilbert Spaces. 2001 . Now, let us fix k^ where r= 0, 1, ... 
- l,inT^sothatKerp= {k, :r = 0, 1,... 1}. 
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